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Within the framework of geodetic brane gravity, the Universe is described as a 4-dimensional 
extended object evolving geodetically in a higher dimensional fiat background. In this paper, by 
introducing a new pair of canonical fields {A, Pa}, we derive the quadratic Hamiltonian for such a 
brane Universe; the inclusion of matter then resembles minimal coupling. Second class constraints 
enter the game, invoking the Dirac bracket formalism. The algebra of the first class constraints 
is calculated, and the BRST generator of the brane Universe turns out to be rank-1. At the 
quantum level, the road is open for canonical and/or functional integral quantization. The main 
advantages of geodetic brane gravity are: (i) It introduces an intrinsic, geometrically originated, 
, 'dark matter' component, (ii) It offers, owing to the Lorentzian bulk time coordinate, a novel solution 

■ to the 'problem of time', and (iii) It enables calculation of meaningful probabilities within quantum 

cosmology without any auxiliary scalar field. Intriguingly, the general relativity limit is associated 
£NJ , with A being a vanishing (degenerate) eigenvalue. 
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I. INTRODUCTION 

(N 

^ I . Geodetic Brane Gravity (GBG) treats the universe as an extended object (brane) evolving geodetically in some flat 
' background. This idea has been proposed more than twenty years ago by Regge and Teitelboim ('General Relativity 
a la String') Q , with the motivation that the first principles which govern the evolution of the entire universe cannot 
be too different from those which determine the world-line behavior of a point particle or the world-sheet behavior of 
' a string. 

Geometrically speaking, the 4-dimensional curved space-time is a hypersurface embedded within a higher dimen- 
sional flat manifold. Following the isometric embedding theorems Q, at most N = |n(n + 1) background flat 
dimensions are required to locally embed a general n-metric. In particular, for n = 4, one needs at most a 10 
dimensional flat background. This number can be reduced, however, if the n-metric admits some Killing-vector fields. 
$H . In Regge and Teitelboim (RT) model, the external manifold (the bulk) is flat and empty, it contains neither a 
gravitational field nor matter fields. Other models were suggested, where the external manifold is more complicated 
, p9|"p3t, it may be curved and contain bulk fields which may interact with the brane. RT action, therefore, does not 
contain bulk integrals, it is only an integral over the brane manifold, which may include the scalar curvature (7?™), a 
rS | constant (A), and some matter Lagrangian (C mat ter) FL 



S = J \m^ 7zn + A + Cmatter ) ~ lxdT (1) 

The Geodetic Brane has two parents: 

1. General relativity gave the Einstein-Hilbert action, which makes the geodetic brane a gravitational theory. 

2. Particle/ String theory gave the embedding coordinates y A (x) ^as canonical fields, and this will lead to geodetic 
evolution. The 4 dimensional metric is not a canonical field, it is just being induced by the embedding g^ y {x) = 
ilABy A l _ l {x)y B ^(x). 



1 The n = 1 brane is a particle, it has 1Z 1 = and A is the mass of the particle. The n = 2 brane is a string, it's curvature 
1Z 2 is just a topological term, and A is the string tension. The brane universe n = 4 includes both the scalar curvature 1Z 4 , 
and the cosmological constant A. 

2 We denote the embedding space indices with upper-case Latin letters, spacetime indices with Greek letters, and space indices 
with lower-case Latin letters, tjab is the Minkowski metric of the embedding space. 
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Due to the fact that the Lagrangian ([!]) does not depend explicitly on y A , but solely on the derivatives through the 
metric, the geodetic brane equations of motion are actually a set of conservation laws 



{TV"' - \^ v Tl - %-kGT» v ) V \ 



= 0. 



(2) 



Eq. (||) splits into two parts, the first is proportional to y A ^ and the second to y A ^ v ■ Since the 4-dimensional covariant 
derivative of the metric vanishes g^-x = 0, one faces the embedding identity tjabV -xV-fiu = 0- Therefore, the first 
and second covariant derivatives of y , viewed as vectors in the external manifold, are orthogonal, and each part of 
Eq.(||) should vanish separately. The part proportional to y A ^ implies that T?^ = 0. The second part is the geodetic 
brane equation 







(3) 



• The matter fields equations remain intact, since the matter Lagrangian depends only on the metric. 

• Energy momentum is conserved. This is a crucial result, especially when the Einstein equations are not at our 
disposal. 

• Clearly, every solution of Einstein equations is automatically a solution of the corresponding geodetic brane 
equations. But the geodetic brane equations allow for different solutions J2|. A general solution of eq.(||) may 
look like 



TV LV - -g^TZ - 8nGT^ 
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fJLV 
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(4a) 
(4b) 



The non vanishing right hand side of eq.(^a|) will be interpreted by an Einstein Physicist as additional matter, 
and since it is not the ordinary it may labeled Dark Matter pqj. 

It has been speculated, relying on the structural similarity to string/membrane theory, that quantum geodetic brane 
gravity may be a somewhat easier task to achieve than quantum general relativity (GR). The trouble is, however, 
that the parent Regge-Teitelboim Hamiltonian has never been derived! 

In this paper, by adding a new non-dynamical canonical field A we derive the quadratic Hamiltonian density of a 
gravitating brane-universe 



H = N' 



y\k 



N 



8nG 
2\/h 



8^G 



(A + -ft (3) ) + J P9(* - XI^QP 



(5) 



The derivation of the Geodetic Brane Hamiltonian is done here in a pedagogical way. In section [n] we translate 
the relevant geometric objects to the language of embedding. Each object is characterized by its tensorial properties 
with respect to both the embedding manifold and the brane manifold. We embed the ADM formalism || in a 
higher dimensional Minkowski background, the 4 dimensional spacetime manifold (V4)is artificially separated into a 
3-dimensional space-like manifold (V3) and a time direction characterized by the time- like unit vector orthogonal to 
V3. For simplicity we restrict ourselves to 3-dimensional space-like manifolds with no boundary (either compact or 
infinite), while the appropriate surface terms should be added when boundaries are present [||. 

Section [III| is the main part of this Paper, where we derive the Hamiltonian. We first look at an empty universe 
with no matter fields, we present the gravitational Lagrangian density as a functional of the embedding vector y A {x), 
and derive the conjugate momenta Pa{x). Reparametrization invariance causes the canonical Hamiltonian to vanish, 
(in a similar way to the ADM-Hamiltonian and string theory), and the total Hamiltonian is a sum of constraints. We 
introduce a new pair of canonical fields A, Pa and make the Hamiltonian quadratic in the momenta. Following Dirac's 
procedure H we separate the constraints into 4 first-class constraints (reflecting reparametrization invariance), and 2 



3 The geodetic factor yf^, 



^^3cV B ^y C v replaces yf^ in case the embedding metric is not Minkowski. 
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second-class constraints (caused by the 2 extra fields). We define the Dirac-Brackets and eliminate the second-class 
constraints. The final algebra of the constraints takes the familiar form of a relativistic theory, such as: The relativistic 
particle, string or membrane. 



In section 



. [V we discuss the inclusion of arbitrary matter fields confined to the four dimensional brane. The algebra 
of the constraints remains unchanged, while the Hamiltonian is simply the sum of the gravitational Hamiltonian and 
the matter Hamiltonian. 

In section 0thc necessary conditions for classical Einstein gravity are formulated, they are 

• A must vanish. 

• The total (bulk) momentum of the brane vanishes. 

Section [v| deals with quantization schemes. We can use canonical quantization by setting the Dirac Brackets 
to be commutators {, } D — ► ]. The wave-functional of a brane-like Universe [ j^J is subject to a Virasoro-type 
momentum constraint equation followed by a Wheeler-de Witt-like equation (first class constraints), the operators are 
not free, but are constrained by the second class constraints as operator identities. Another quantization scheme is 
the functional integral formalism, where we use the BFV |ll| formulation. The BRST generator |L7| is calculated, 
and the theory turns out to be rank 1. This resembles ordinary gravity and string theory as oppose to membrane 
theory, w here the rank is the dimension of the underlying space manifold. 



Section VII Geodetic Brane Quantum Cosmology is demonstrated. We apply the path integral quantization to the 
homogeneous and isotropic geodetic brane, within the minisuperspace model. A possible solution to the problem of 
time arises when one notices that while in GR the only dynamical degree of freedom is the scale factor of the universe, 
GBQC offers one extra dynamical degree of freedom (the bulk time) that may serve as time coordinate. 

Definitions, notations and some lengthy calculations were removed from the main stream of this work and were put 
in the appendix section. 



II. THE GEOMETRY OF EMBEDDING 



In this section we will formulate the relevant geometrical objects of the V4 and V3 manifolds in the language of 
embedding. Let our starting point be a flat m-dimensional manifold M. , with the corresponding line-element being 



ds 2 = r/ABdy A dy B . (6) 

Hypersurfaces : An embedding function y A {x'- i ) (/i = 0,1,2,3) defines the 4 dimensional hypersurface V4 pa- 
rameterized by the 4 coordinates x^ . The V4 tangent space is spanned by the vectors y A . (The V3 hypersurface 
and tangent space are defined in a similar way). The induced 4-dimensional metric is the projection of tiab 
onto the V4 manifold: — T)ab% I uV B v - Choosing a time direction t and space coordinates x % (i = 1,2,3), the 
induced 4-dimensional line-element takes the form 

ds 2 = VABiyidx* + i/'dhd/'Ur* + y B dt) , (7) 

The various projections of the metric tjab onto the space and time directions are denoted as the 3-metric hij, 
the shift vector Ni, and the lapse function N 

VAbV^Vj = hij , (8a) 
VABviy 3 = Ni , (8b) 
VABV A y B = NiN* - N 2 . (8c) 

These are not independent fields (as in Einstein's gravity), but are functions of the embedding vector y A . 
Nevertheless, it is a matter of convenience to write down the induced 4-dimcnsional line-element in the familiar 
Arnowitt-Deser-Misner || (ADM) form 

ds 2 = -N 2 dt 2 + h ij {dx i + N l dt)(dx ] + N j dt) . (9) 

The vectors (y A ,y A ) span the 4-dimensional tangent-space of the V4 spacetime manifold, while y A span the 
3-dimensional tangent-space of the V3 manifold. Using h? 3 as the inverse of the 3-metric h'^hjk = 6\, one can 
introduce projections orthogonal to the V3 manifold with the operator 
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©s = Si - y A a h ab y B ,b . (10a) 

©£©£ = &£■ . (iob) 

Now, any vector v A can be separated into the projections tangent and orthogonal to the V3 space 

v A = «f + «f = ^2/b,6/i o6 2/^ + « B (^ - • (11) 

An important role is played by the time-like unit vector orthogonal to V3-space yet tangent to V4-spacetime, 

n A ^±(y A ~N*y A ) = ±y B e A , (12a) 

VABV A t n B = , (12b) 
j]ABn A n B = -1 . (12c) 

The tangent space of the embedding manifold M. is spanned by the vectors: y A , n A and L A (i = 1,2, 3 p = 
1, .., m — 4). The vectors L A are chosen to be orthogonal to y A , n A and to each other. 

• Curvature : The connections on the underlying V3 are = r]ABy A jy B h kl , this way, the covariant derivative of 
the 3-metric vanishes /lyife = (the stroke denotes 3-dimensional covariant derivative). As a result, one faces 



the powerful embedding identity 

~"i\ijV,k 



VA B yty B k = . (13) 



The vectors yy- . are orthogonal to the V3 tangent space and may be written as a combination of n and L A 



y%=n A K l]+ L A % . (14) 

. i.i A 

K '^-^( N *\i + N i\*-^) = - l i^<i ,,U ■ [ir>) 



The projection of yf- ■ in the n A direction is the extrinsic curvature of the V3 hypersurface embedded in V4 



' ' • ■ " 'H'ij \ ,. ,,.\ ,,u 

-A 



The coefficient fi« is the extrinsic curvature of V3 with respect to the corresponding normal vector . 

The intrinsic curvature of the V3 manifold is also related to the second derivative of the embedding functions 
y A j- The 3-dimensional Riemann tensor is 

4% = Vab G4»f H - yfayfji )■ (16) 

For convenience we define the y^-independent symmetric tensor 

V AB = (h^h ab - K a h? h )y A 3 y B ah . (17) 

Checking the indices, ^ AB is a tensor in the embedding manifold, but a scalar in V3 space. The trace of ^ A B is 
simply the 3-dimensional Ricci scalar IZ^ — t/ab ^ AB ■ Looking at eq.([l5|), one can easily check that 

Hy* = , (is) 



and \& as an operator has at least 3 eigenvectors with vanishing eigenvalue. Using definitions (|17J,y_5|), the 
contraction of ^ twice with n A is related to the extrinsic curvature 

K\K] - KijKV = y AB n A n B = ^ AB y A y B . (19) 
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III. DERIVING THE HAMILTONIAN 



The gravitational Lagrangian density is the standard one 

1 



16nG 



-9 



Up to a surface term, it can be written in the form 

1 



C = 



16ttG 



NVh K i3) - {K\K\ - KijK ij ) 



(20) 



(21) 



Here, denotes the 3-dimensional Ricci scalar, constructed by means of the 3- metric (^a|), whereas Kij (|l5| ) is 
the extrinsic curvature of V3 embedded in V4. Using the tensor ^i AB ( |l7j ) one can put the Lagrangian density ( |2lD in 
the form 



C = 



16ttG 



(22) 



As one can see, the Lagrangian (22) does not involve mixed derivative y A or second time derivative y . The first 
derivative y A appears either explicitly or within N . Therefore the Lagrangian 



£{y,y,y\i,y\ij) 



(23) 



is ripe for the Hamiltonian formalism. 

The momenta Pa conjugate to y A is simply 



Pa{x) 



5L 



Sy A (x) ~ 16nG 



+ ^scy B y c 



ON 
dy A 



N 



^ABV 



(24) 



Using eq.(8_b,3cj) to get = -n^, while eq. (|l8|) tells us that j^^abV 3 — ^ABn B , the momentum ( gjj ) becomes 



P A = — 



n c 



i A + 2»£ A n B 



} 



(25) 



The next step should be :" Solve eq.(^5|) for y A (y, P, y\j, y\jj) " ■ But cq.(p5|) involves only n A , so one would like to solve 
eq.(|2^) for n A (P, y, yu, yuj) first, and then to solve eq.(12a) for y A 



y A = Nn A + N'lfi 



(26) 

This looks innocent but even if one is able to solve eq.(p5|) for n A , any attempt to solve eq. (pb|) for N % {n, y, yu) and 
eq.(p^) for N(n, y, y^) will lead to a cyclic redefinition of N l and N. This situation is similar to other reparametrization 
invariant theories (such as the relativistic free particle, string theory etc.) and simply means that we have here 4 x V3 
primary constraints 



i] AB n A n B + 1 = 
VAByin B = . 



(27a) 
(27b) 



The constraints should be written in terms of canonical fields (jj a ,Pa)- So one should solve eq.(^5|) for n A (P), and 
then substitute in the above constraints. Any naive attempt to solve eq.(p5|) for n A (y,P) falls short. The cubic 
equation involved rarely admits simple solutions. To 'linearize' the problem we define a new quantity A, such that 



P- 



^ {*-\I) B n B 



8ttG 



(28) 



• Comparing eq.(|25|) with eq.(|28j), the definition of A is actually another constraint 

n A ^ A B n B +K {3) +2A = . 



(29) 
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• An independent A comes along with its conjugate momentum P\. A is not a dynamical field therefore one faces 
another constraint 



Pa =0 



(30) 



Assuming A is not an eigenvalue of , we solve <M) for n A (y/h, P, A) and find 



8ttG 

~7T 



(* - xiy 



p l 



(31) 



At this point we have 6 x V3 primary constraints ( 27a , 27h . 29| , |30j ) . We will follow Dirac's way J?! to treat the constrained 
field theory we have in hand. 

First we will write down the various constraints in term of the canonical fields (y (x), Pa{x), X(x), P\(x)) : 



8ttG 
2Vh 



' \[h ' 
8^G 



(A + TZ (3) ) + P0(tf - XI^QP 



P A ~0 , 



8ttG 
2Vh 



8^G 



pe(* - xiy 2 &p 



(32a) 

(32b) 
(32c) 

(32d) 



Notations 



• We use shorthanded notation to simplify the detailed expressions, F ■ G = F a Ga where F and G are vectors in 
the embedding space, and P(W - AP)~ 2 P = Pa [(* - XI)- 2 ] AB P B 

• We adopt Dirac's notation <f> pz for weakly vanishing terms. 

• The embedding functions y (x) and A(a:) are scalars in the V3 manifold. Their conjugate momenta Pa(x), P\{x) 
are scalar densities of weight 1. For convenience we normalize all constraints to be scalars in the embedding 
space, and scalar /vector densities of weight 1 in V3. This way, the Lagrange multipliers are of weight 0. 



(33) 



is based on the constraint ( |27b|) but it takes into account the embedding identity (jig) 

<Pk = y\k- P = _ ^~^y|fe(* ~ A/ ) n = ^~^y\ k ' n ~ • 



8vrG" 



• 05 is based on the constraint (27a), but we added the projection operator 9 (10b) in front of P. This step 
simplifies the final algebra of the constraints, and brings it to the f amil iar form of a relativistic theory. Inserting 
O in front of P, is equivalent to adding terms proportional to <j>k (|32b|), since 



6£P* = (8* - v%K»v B , h )P» = P A - y A a h ab fo . 



(34) 



6 is also a combination of the constraints (Eg), (27b) and (27a), chosen such that 



dX 







(35) 



• See appendix ^ for the definitions of functional derivatives and Poisson brackets. 
In a similar way to other parameterized theories, the canonical Hamiltonian density vanishes 

n c = y A P A - C w . 
This means that the total Hamiltonian is a sum of constraints 



(36) 



G 



H = / d A xu m {x)ct> m {x) 



(37) 



The constraints ( |32| ) should vanish for all times, therefore their PB with the Hamiltonian should vanish (at least 
weakly). This imposes a set of consistency conditions for the functions u m (x) 

<j> n {x) = {<f> n (x),H} = {(f> n (x), J d 3 zu m {z)<j) m {z) 

d 3 zu m (z){ ( j >n (x),cl )m (z)}^0 . (38) 
The basic Poisson brackets between the constraints are calculated in appendix pL and in general has the form 



{,}« 


4>o{z) 




4> 4 {z) 


05 W 


<f> (x) 











^ 


<t>k{x) 











&5,\\k<>(x - z) 


04 (a;) 











~<l>5,\S(x - z) 






-05,AA|A;(5(a; - z) 




[F l (x)+F^(z)}S ll (x^z) 



(39) 



The exact expressions for a and F % appears in appendix [b|. Now, insert the PB between the constraints (|3^) into the 
consistency conditions (pq) to determine u m (x) 



dX 



{x)u 5 {x) w =*> u 5 (a;) = , 



{00 (a), -ff} ~ =>• w°(a;) = -N(x) arbitrary , 
{0 fe (a;),il} « =>• u k {x) = N k (x) arbitrary , 
"8vrG 



Mx),H} 



dh 



2\fh 



(z, x)N(z) - <p 5 , x X\ k S(x - z)N k (z) + <p 5 , x S(x - z)u 4 (z) 



u 4 (x) = N k X {k (x) - (f>l\(x) j d 3 z ^-j=a(z,x)N(z) 



8ttG 



(40) 

(41) 
(42) 



(43) 



The first class Hamiltonian is then 



H = J d 3 x[N k [y ]k ■ P + X ]k P x ] 



N 



8ttG 



87TG 



2Vh 

+ I d 3 za{x,z)4)-^(z)Px{z) 



(44) 

As one can see, at this stage we have in the Hamiltonian 4 arbitrary functions N,N k (Lagrange multipliers). This 
means we have 4 first class constraints reflecting the reparametrization invariance (4-dimensional general coordinate 
transformation) 



•■Po 



8ttG 
2\/h 



{^) 2 {X + K {3) ) + P6(* - A7)- x eP 



8ttG 

+ / d 3 za(x,z)<j)~ 1 x (z)Px(z) 
fk =y\ k P + X\ k P x w . 



(45a) 
(45b) 



And we are left with 2 second class constraints, reflecting the fact that we expanded our phase space with two extra 
fields A and Pa 



h = 04 = Pa « 
8vrG 



2Vh 



8^G 



pe(* - xiy 2 @p 



(46a) 
(46b) 
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Using the classical equation of motion for y A (x), 

y A (x) = {y A (x),H}^N k y lk 



N^(*-XI)^P, 



(47) 



one can identify the lapse function ( |8q) and the shift vector (8b) with N,N k respectively. Thus, recover the nature 
of the lapse function and the shift vector as Lagrange multipliers only at the stage of the solution to the equation of 
motion, not as an a priori definition. 

We would like to continue along Dirac's path B, and use Dirac Brackets (DB) instead of Poisson Brackets (PB). 
The DB are designed in a way such that the DB of a first class constraint with anything is weakly the same as the 
corresponding PB, while the DB of a second class constraint with anything vanish identically. Using DB, we actually 
eliminate the second class constraints (the extra degrees of freedom). The DB are defined as 



d 3 x / d 3 z{A,e m (x)} P C- 1 n (x,z){e n (z),B} 1 



{A,B} D ^{A,B} P 

Where (a;,z) is the inverse of the second class constraints PB matrix 

C mn (x,z) = {9 m (x),9 n (z)} . 
In our case, C mn (x, z) is simply the 2x2 bottom right corner of ( |39| ) 

/ 



(48) 












[F i (x) + F i (z)] 6 {i (x-z) 



1,2 



(49) 



When dealing with field theory, the matrix C mn is generally a differential operator, and the inverse matrix is not 
unique unless one specifies the boundary conditions. We choose 'no-boundary' as our boundary condition, therefore, 
integration by parts can be done freely, and the inverse matrix is 



The resulting DB are 

{A,B} D ={A,B} P 



d 3 x 



((tr 2 n*) + < 




S\i(x - z) 




KdXj 


(x)S(x - z) ^ 




ydxj 


-1 

(x)S(x — z) 


) 



(50) 



SA SB SA , 



SB 



SX(x) 5X(x) 



dX 



(x) 



SA 

SX(x] 



{^{x),B} + {A^ b {x)} 



5X{x) ni 5X{x) 
SB 



SX(x) 



(51) 



This way, from now on, one should work with DB instead of PB and take the second class constraints to vanish 
stron gly. This will omit the parts proportional to P\ from the first class constraints ( 454 45b) and recover the original 
form ( |32a|j32b| ). 

The algebra of the first class constraints takes the familiar form Q| of a relativistic theory 



{(f>o(x),(f> {z)} D = [h ij (j)i{x) + /i y i (z)]%(2;- z) 
{(f>o(x),(f> k (z)} D = (j> (z)5\ k (x - z) 
{(f>k(x),<f>i(z)} D = <pi(x)5\ k (x - z) + (f> k (z)S\i(x - z) 

The final first class Hamiltonian of a bubble universe is 



H = J d 3 x^N k y lk -P-N 



8ttG 
2y/h 



'VV 

8nG 



(A + ^ (3) ) + P9(* - A7) _1 6P 



(52a) 
(52b) 
(52c) 



(53) 



At this stage, we have a first class Hamiltonian composed of four first class constraints, and accompanied with 
two second class constraints. The algebra of the first class constraints is the familiar algebra of other relativistic 
theories. Before moving on to quantization schemes we would like to study two more classical aspects: what happens 
if the action includes brane matter fields, and what is the relation between Einstein's solutions to the geodetic brane 
solutions. 
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IV. INCLUSION OF MATTER 



The inclusion of matter is done by adding the action of the matter fields to the gravitational action 

1 



S= d 



J4, 



16ttG 



(54) 



The matter Lagrangian density depends in general on some matter fields, but also on the 4-dimensional metric g^ v . 
The dynamics of the matter fields is actually not affected by the exchange of the canonical fields from to y A , and 
one expects the same equations of motion or the same 'matter' Hamiltonian density. On the other hand the momenta 
Pa gets a contribution from the matter Lagrangian 

AP A = ^j^ = Vh[T nn n A -^T ni y A i \ . (55) 

This contribution depends on the various projections of the energy-momentum tensor 

rpfiu _ 2 SC ma tter (561 

T nn is the matter energy density, or the projection of the energy-momentum tensor twice onto the n A direction 
T nn = {T^ v y A y^j uatib- While in T n i the energy-momentum tensor is projected once onto the n A direction and 
once onto the V3 tangent space. T n i = (T^y A ^y B ^ nAl)B,i- See appendix |c| for some examples of matter Lagrangians, 
Hamiltonians and the corresponding energy-momentum tensor projections. 
The momenta Pa (p5| ) is now changed to 

P A = { [^ (3) + n B * BC n c - 167rGT n „] n A + 2* A n B + 16nGT ni h^y A } (57) 

Following the same logic that lead us from Eq. ( p5| ) to the introduction of A (^9|) , we will define A as 

n A ^ A B n B + K {3) - 16nGT nn + 2A = . (58) 

The effects of matter are thus A — > A + 87rGT„„, P a — > P A — ^/hT ni h % ^y A but OP is unchanged. The constraints are 
modified as follows 

4>Q > 00 - VhT nn (59) 

(/>k — > <i>k + VhT nk (60) 

Thus the Hamiltonian is changed to 

Hg — * Hg + j d 3 xVh,[N k T nk + NT nn ] =H G + H m , (61) 

Where H m is the matter Hamiltonian, calculated in terms of the matter fields alone as shown in appendix |^. The 
algebra of the constraints ( js2| ) remains unchanged under the inclusion of matter, where the PB now include the 
derivatives with respect to matter fields as well. 

V. THE EINSTEIN LIMIT 

In some manner Regge-Teitelboim gravity is a generalization of Einstein gravity. Any solution to Einstein equations 
is also a solution to RT equations (||). We will derive here the necessary conditions for a RT-solution to be an Einstein- 
solution. 

• First, we use a purely geometric relation 

2G nn =ft( 3) +n B ^ BC n c , (62) 
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where G nn is the Einstein tensor twice projected onto the n^-direction. The constraint associated with the 
introduction of A (|5|) is 



- 2A = Tl^ + n B *Bcn c - WnGT nn = 2{G nn - 8irGT nn ) 
The Einstein solution of the equation is therefore associated with 



A = 



(63) 



(64) 



As was shown in Eq.(|lJ), \1/ has a degenerate vanishing eigenvalue. Therefore Einstein case with A = 0, will not 
allow for the essential (\& — A/) -1 . One can not impose A = as an additional constraint (as was proposed by 
RT [0), but only look at it as a limiting case. 

• Second, we use the projection of the Einstein tensor once onto the n^-direction and once onto the V3 tangent 
space G n i 



G,„l>''!l\ = ~^ A B n B - (v A j [Kh ij - , 



in eq.(^) and put the momentum P A in the form 
Vh 



P A = 



8nG 



(G nn - 8irGT nn )n A ~ (G m ~ 8nGT ni )h ij y A + ( y l , •: A /' ' *' - K**)) 



(65) 



(66) 



It is clear that if Einstein equations G nn — 8nGT nn and G n i = 8nGT n i are both satisfied, the momentum Pa 
makes a total derivative such that 



d 3 xP A = 



(67) 



The total momentum § cPxPa is a conserved Noether charge since the original Lagrangian does not depend 
explicitly on y A 



■ A - * d 3 xP A = const. 



(68) 



The universe, as an extended object, is characterized by the total momentum /i A . The necessary condition for 
an Einstein-solution is a vanishing /i A . 



d 3 xP A = 



(69) 



• The condition (^) simply tells us that the total 'bulk' momentum of the universe vanishes. This motivates us 
to use a new coordinate system for the embedding, namely, the 'center of mass frame' + 'relative coordinates'. 
As relative coordinates we will use the derivatives y A . This has a direct relation to the metric and therefore, we 
expect the equation of motion to resemble Einstein's equations. The new system and the calculations appear in 
appendix [d| 



VI. QUANTIZATION 

The treatment so far was classical, but the derivation of the Hamiltonian and the construction of the various 
constraints are the ingredients one needs for quantization. In the following sections we will describe two quantization 
schemes, canonical quantization and functional integral quantization. 
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A. Canonical Quantization 



Dirac's procedure leads us towards the canonical quantization of our constrained system. The following recipe was 
constructed by Dirac H] for quantizing a constrained system within the Schrodinger picture 

• Represent the system with a state vector (wave functional) . 

• Replace all observables with operators. 

• Replace DB with commutators, {, } D — > ih [, ] 

• First class constraints annihilate the state vector. 

• Second class constraints represent operator identities. 

• Since the commutator is ill defined for fields at the same space point, one must place all momenta to the right 
of the constraint. 

• First class constraints must commute with each other. This ensures consistency, and may call for operator 
ordering within the constraint. 

In our case, we can use the coordinate representation. The state vector is represented by a wave functional $[y]. The 
DB (commutator) between y A and Pb is canonical, therefore, these operators can be represented in a canonical way 

y A {x)^y A {x) 



Pa(x) 



-ih- 



6y A (x) 



The operator P\ vanishes identically. The DB of A with y ,Pb are not canonical, therefore the operator A must be 



expressed as a function of y A ,Ps ■ This can be done with the aid of the second class constraint ( 46b ). 

The first class constraints as operators must annihilate the wave functional. These constraints are recognized as 



1. The momentum constraint (J45b|) 



a 5$ n 



(70) 



which simply means that the wave functional is a V3 scalar and does not change its value under reparametrization 
of the space coordinates. This can be shown if one takes an infinitesimal coordinate transformation 



c fe + e k 



y A {x)-^y A {x)+e*y%{x) , 



The wave functional is unchanged if and only if the momentum constraint holds. 

2. The other constraint is the Hamiltonian constraint, and up to order ambiguities, the equation is the analog to 
the Wheeler de-Witt equation 



8ttG 
2Vh 



8^G 



(\ + K^)(x)-h 2 ((^-A/)- 1 ] (x) 



5y A (x)6y B (x) 



$[y] = 



(71) 



It is accompanied however, with the operator identity 

2 



8ttG 
2Vh 



8^G 



pe(#- \iy 2 op 



= 



(72) 
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B. Functional Integral Quantization 



Calculating functional integrals for a constrained system is not new. This was done for first class constraints by 
BFV |lf|, And was generalized for second class constraints by Fradkin and Fradkina flG|| . 

The first step is actually a classical calculation, that is, calculating the BRST generator |L7|. For this calculation 
we will adopt the following notations: 

• The set of canonical fields will include the Lagrange multipliers N 1 - 1 = (N,N l ) that is Q A — (y A , A, A"') T , 
and the corresponding conjugate momenta ILa = {Pa, P\, 7!"^). The Lagrange multipliers are not dynamical, 
therefore, the conjugate momenta must vanish. This doubles the number of first class constraints G a — (vr^, <j> v ). 

• For each constraint we introduce a pair of fermionic fields i] a — (p^,^) 1 ", and the conjugate momenta V a — 
(c v ,p v ). (In our case, all constraints are bosonic, therefore the ghost fields are fermions). 

• Each index actually represent a discrete index and a continuous index, for example, y A = y (x). The summation 
convention is then generalized to sum over the continuous index as well 

N^= [ d 3 xN»(x)^(x) . (73) 



• We use Dirac Brackets as in ([3lJ), but the Poisson Brackets are generalized to include bosonic and fermionic 
degrees of freedom 

Where (q,p) is the set of canonical fields including the fermionic fields, r, I denote right and left derivatives 

Jn d r R J , d l R 

dR = —dq = dq— . (75) 



And the fermionic index is 



1 if R is a boson 
n R — S 1 if r i s a fermion ' 



Let us now calculate the structure functions of the theory. The first order structure functions are defined by the 
algebra of the constraints {G ai Gb}D — G c U^ b . It is only the original constraints, (not the multipliers momenta), that 
have non vanishing structure functions (|52|). 

{(ZU D = (0 f<PwM*)U*,(x,z,w) ) ' (77) 

and the relevant first order structure functions are 

U^(x, z, w) = [Sffih™ (8(w -x) + 5(w - z)) + 5ffi6(w - z) + 8ffi5{w - x)\ S tk {x - z) (78) 

(Generally, one should also look at {H ,G a }£> = G^V^i but here H — 0). The second order structure functions are 
defined by the Jacobi identity A({{G a , G{,}d, G c }d) = 0, where A means antisymmetrization. Using the first order 
functions (|7|) one gets A(G d [{U* b , G C } D + U* c Ul b }) = 0. This equation is satisfied if and only if the expression in the 
square brackets is again a sum of constraints 

A{{U d ab , G C } D + U d ec Ul b ) = GfUH. (79) 

The second order structure functions U[ b d c are antisymmetric on both sets of indices. In our case, the second order 
structure functions vanish, and the theory is of rank 1. This resembles ordinary gravity and string theory as oppose 
to membrane theory, where the rank is the dimension of the underlying space manifold. The BRST generator of a 
rank 1 theory is given by Q — G a rj a + ^'PcU ° b r] b ri a . Here it is 
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n = / d 3 x[ir^ + <t>^ + h kl p k c° lC ° + p^c^x). (80) 



The main theorem of BFV [jl5[ is that the following functional integral does not depend on the choice of the gauge 
hxing Fermi function W : 

= J VQ A VU A Vr] a VVa M exp[i J dt(U A Q A + VJf - . (81) 

Where M = S{8 1 )6(8 2 ){detC mn ) 1 / 2 is taking care of the second class constraints, and, since the canonical Hamiltonian 
vanishes, = —{ty,Q}D- 

The determinant of C mn for compact space manifolds, is calculated in a simple way in appendix ^. 

VII. AN EXAMPLE: GEODETIC BRANE QUANTUM COSMOLOGY 

In the following example we would like to implement GBG to cosmology, and in particular to quantum cosmology. 
Detailed examples and calculations can be found in |^,^9|, here we will just focus on global characteristics of the 
Feynman propagator for a geodetic brane within the minisuperspace model. Attention will be given to the differences 
between 'Geodetic Brane Quantum Cosmology' and the standard 'Quantum Cosmology'. 

The standard and simple way to describe the cosmological evolution of the universe is to assume that on large scales 
the universe is homogeneous and isotropic. The geometry of such a universe is described by the Friedman-Robertson- 
Walker (FRW) metric 

ds 2 = -N 2 (t)dt 2 + a 2 (t)dtt 2 3 , (82) 

where N(t) is the lapse function, a(t) is the scale factor of the universe, and 

dQl = dip 2 + X 2 {^)dnl (83) 

is the line element of the 3 dimensional spacelike hypersurface which is assumed to be homogeneous and isotropic, dfl 2 
is the usual line element on a 2 sphere, and x(V0 = sin-0 / ip / sinh^ if the 3 space is closed, flat or open respectively. 
In General Relativity, the components of the metric are the dynamical fields, the lapse function N(t) is actually a 
Lagrange multiplier, and the only dynamical variable is the scale factor a(t). This model is called minisuperspace, 
since, the infinite number of degrees of freedom in the metric is reduced to a finite number. The remnant of general 
coordinate transformation invariance, is time reparameterization invariance, that is, the arbitrariness in choosing 
N(t). The usual and most convinient gauge is N = 1. 

In GBG the situation is quite different. First, one has to embed the FRW metric (|82|) in a flat manifold. The 
minimal embedding of a FRW metric calls for one extra dimension. We will work here, for simplicity, with the closed 
universe x = sini/>. The embedding in a flat Minkowski spacetime with the signature (— , +, +, +, +), is given by M 



V A -( T[t) r ) z 1 - 



/ sin ip sin 8 cos < 
sin ip sin 8 sin c, 
sin ip cos 9 

\ cost/ 1 



(84) 



The lapse function is given by N(t) = \/T 2 — a 2 , it is not a Lagrange multiplier, but it depends on the two dynamical 
variables: the scale factor a(t) and the external timelike coordinate T(t). Time reparameterization invariance is, 

naturally, an intrinsic feature of \/T 2 — a 2 dt, but, no gauge fixing is allowed here, since both T(t) and a(t) are 
dynamical. The gravitational Lagrangian (^2|), after integrating over the spatial manifold is 

L = a[3Na- — 

a = is a scaling factor, for convinience we will set a = 1. The key for quantization is of course the Hamiltonian. 
One can derive the Hamiltonian directly from the Langrangian (|S5|), or, to use the ready made Hamiltonian ( |53| ) and 
just insert the 'minimized' expressions for the embedding vector, and the conjugate momenta. 
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A. Minisuperspace Hamiltonian 



The first step is to introduce the coordinates and conjugate momenta. The general embedding vector y A is replaced 
by the dynamical degrees of freedom a(t) and T(t), while the spatial dependence is forced by the expression (|8J). It 
is expected that the conjugate momenta will have two degrees of freedom P a (t), Pr(t), the delicate issue is the spatial 
dependence of the momenta. Our choice is 



Pa 



P T (t) 



sin i/j sin 9 
8ttG 



(86) 



the factor sin 2 ip sin 9 was inserted in order to keep the moment a a 3-dimensional vector density. The spatial de- 
pendence is through z I (x) such that the momentum constraint ( |32b| ) vanishes strongly. And, the normalization is 

/ d 3 xy A pA = c(ctP a + TPt). In addition, we set A 
Inserting these expressions into the constraints ([ 
first class constraint 



sin" ?/; sin 
8^G 



= A(t) and P x = P x (t) 

a) and integrating over spatial coordinates, one is left with one 



+ a 3 A + 



P 2 



P 2 

T 



a 3 A 6a — a 3 A 



aPy 



(87) 



and two second class constraints 



h =Px 



The Dirac brackets (51) arc defined as 



a 3 - 



p2 
T 



a 3 P 2 



i 3 X 2 (6a - a 3 A) 2 



{AB} D ={A,B} } 



and the minisuperspace Hamiltonian is 



H 



-N 
~~2~ 



p2 
T 



a b P 2 



i 6 A 3 (6a~a 3 \y 



6a + a 3 A + 



P 2 

a 3 A 6a 



aP x 



(89) 



(90) 



We would like to focus on the Feynman propagator JTJ] K(af 1 Tf,tf 1 a i ,T i ,ti) for the empty geodetic brane universe. 
Although the empty universe is a non-realistic model for our universe, the calculation of the propagator is simple and 
it demonstrates some of the main features and advantages of Geodetic Brane Quantum Cosmology over the standard 
quantum cosmology models. This propagator is the probability amplitude that the universe is in (a/, Tj) at time tf, 
and it was in (ai,Ti) at time t%. We will use a modified version of BFV integral offered by Senjanovic [pL8[j , where the 
ghosts and multipliers were integrated out. 



K{ a fiTf,tf,ai,Ti,ti) 



dfj, exp 



2iri / dt{aP a +TP T + XP X ) 



dfx = dadP a dTdP T dXdP x 6(<p) S( X ) \{x, <f}\ <*(0i) S(6 2 ) | det({0 m , 9 n })\ 1 ' 2 



(91) 



This propagator is calculated in phase space, where the measure is the Liuville measure dx dp. In addition, the 
measure d/i enforces the constraints (first and second class) by delta functions, it includes an arbitrary gauge fixing 
function the determinants of the Poisson brackets between first class constraints and the gauge fixing function and 
the determinants of the Poisson brackets between second class constraints. Attention should be given to the following 



• The canonical Hamiltonian vanishes, therefore it is absent in the action. 

• The boundary conditions for the propagator determine the values of af,Tf,a,i,Ti, but not the value of A nor 
the values of the momenta. Therefore, the momenta and A must be integrated over at the initial point. 

• The gauge fixing function Xj although arbitrary, must be chosen such that it does not violate the boundary 
conditions nor the constraints. In addition, the Poisson brackets {x, tp} must not vanish. 
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The determinant of the second class constraints Poisson brackets is simply 



\&et{{0 m ,e n })\ 1 ' 2 



d6 2 




Pt a 6 P a 2 


d\ 




a 3 A 3 ' (6a — a 3 A) 3 



(92) 



• Our convention here is a = 1 and Planck constant h — 1 {% = ^— )• 

• In cases where matter is included, the inclusion of matter will affect in a few places. The action will include 
terms like </>7r, an integration over matter fields and momenta will be added, and the first class constraint will 
have a contribution which is simply the matter Hamiltonian tp — > ip + H m (a, <fi, it). All other constraints remain 
intact. 

The calculation of the propag ator (|l]) is carried out in a simple way following Halliwell , and the final propagator 
takes the form 

K±{af,T f ;ai,Ti) = J doj exp \2mu(Tf — Ti) =p 2mw 2 (F(xf) — F(xi))] (93) 

The index of K± and the =F in the exponent refers to the expanding/contracting scale factor, ui is the conserved 
bulk energy (the momentum conjugate to the bulk time coordinate T). Since the value of uj is not fixed at the initial 
condition, one must integrate over lo. One should notice according to eq. (|69j) that the Einstein solution is assiciated 
with u> = 0. The function F(x) is given by 

FY \ _ / 12 [SArcsinx + y/l - x 2 (Ax 5 + 2x 3 - 3a;)] \x\ < 1 , > 

[X) ~ \ sgn(»§ - ± [3sgn(x)Arccosh|x| - Vx^l(4x 5 + 2x 3 - 3x)] 1 < \x\ [ ' 

Where x = (^j) 1/3 . 

Let us now examine the properties of the propagator (|93|). Actually, the propagator is independent of the internal 
time parameter t (a common character of all parameterized theories) , and depends exclusively on the value of a and 
T at the boundaries. 

• The most basic characteristic of a propagator is the possibility to propagate from an initial state to a 
final state through an intermidiate state. For example, the propagator for a non-relativistic particle is 
K(x3, t^; x\, t\) — f dx2K(x3,t3;x2,t2)K(x2,t2',xi,ti). At the intemidiate time t<ii one must integrate over 
X2- It is clear that there is no integration over t^-, t is the evolution parameter, it must be mono- 
tonic £3 > t2 > ti, and integration over t% makes no sense. Another characteristic of the propagator is 
lim t2 ^ tl K(x2, t2] x\, ti) = S(x2 — x\). The situation with parameterized theories is quite different. The propa- 
gator is independent of the internal time, and integration over all dynamical variable diverges. The solution is, 
usually, to let one of the dynamical variables as 'time', and integrate only over the other variables. 

The question is: How does the propagator ( |93| ) behaves at the intermidiate point ? What is the relevant evolution 
parameter and what integrations should be made ? One can check that if a is taken to be the monotonic evolution 
parameter and integration over T at the intermidiate point is done, then the propagator (B3) is well behaved. 



A-(aa, TsjolTx) - / dT 2 J dcj e 2™KT 3 -T 2 )-^(nx 3 )-F(* 2 ))] / ^^[o^-t^-^f^-f^))] 

duo e 2 ^(n-T 1 )^ 2 (F(x 3 )~F(x 1 ))] ^ 

A-(oi,T 2 ;oi,Ti) = / dcue 2 ™^-^ = 6(T 2 - Ti) (96) 



This cannot be done within the standard quantum cosmology models, since there, the only dynamical variable 
is the scale factor a. Such a propagator of only one variable contains no information, it can tell that the varible 
is monotonic. The common solution in standard quantum cosmology is to add another dynamical variable such 
as a scalar field and to use one of them as the evolution parameter. Here we see one of the main advantages of 
Geodetic Brane Quantum Cosmology over the standard models, the problem of time has an intrinsic solution 
as we have one extra degree of freedom which serves as 'time'. 
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The most general wave function that can be generated using the propagator (|93|) is 
*(a,r)= fdwe 2niuT \A(u)e- 2 ™ 2F W +B(u)e 2niu2F W~ . 



(97) 



One can verify that the wave function ( |97| ) (and the propagator (93)) satisfy the corresponding WDW equation 

-t(*)iz ( TTZT^T ) + ^)^F5 *(a,r)=0 (98) 



9a VC(a;)9a; 5 1 ' dT 2 ^ 

Where £(x) = (1 + 2x 2 )\/l — x 2 , and x = (3a(— ift^) -1 ) 1 ^ 3 . Putting —ih-^p — lu and neglecting the term 
proportional to the first derivative ^ , the equation ( |98| ) looks like a zero energy Schrodinger equation 



with the potential 

V u (a) 



r) 2 

-.2 ° 

da 2 



-n z — + V^(a) 



* u (o) = , 



3a s 









^2/3" 
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FIG. 1. The potential V u (a) 



1.4 



The classical turning point is a = cj/3, and the empty brane universe can not expand classically bycond this 
point. The empty universe model is non-realistic, a more realistic model may include some matter fields, or at 
least a cosmological constant. Analysis of the cosmological constant universe can be found in |27}| . 

In accordance with section |v|, one of the necessary conditions for an Einstein solution is Eq. ( |o7|) / d 3 x Pa = 0. 
Within our minisuperspace model, integrating the momenta (^6|) over the spatial manifold one gets J d 3 x Pa = 
(Pt, 0, 0, 0, 0) T , thus the Einstein case is associated with w = 0, and the only classcal regime is a = 0. 
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• The still open question is that of the boundary conditions. In particular \&(a = 0, T) and \&(a — > oo,T). One 
possibility is: "J vanishes at the big bang (a = 0) and W is bounded at a — > oo. This will lead to a; quantization 
cj^ = 8h(n + 1/4) where n is a positive integer. Clearly, the Einstein case oj = is excluded by such quatization 
condition. 



1. In the present model of Geodetic Brane Gravity, the 4 dimensional universe floats as an extended object within 
a flat m dimensional manifold. It can be generalized however, to include fields in the surrounding manifold 
(bulk), this is done by adding the bulk action integral to the action of the brane. The brane will feel those 
bulk fields as forces influencing its motion [122] . The bulk fields may include matter fields or the bulk gravity 



2. In this paper we have derived the quadratic Hamiltonian of a brane universe. The Hamiltonian is a sum of 
4 first class constraints, while 2 additional second class constraints are present. We used Dirac Brackets and 
found the algebra of first class constraints to be the familiar one from other relativistic theories (such as string, 
membrane or general relativity). The BRST generator turns out to be of rank 1. 

3. Geodetic brane gravity modifies general relativity, and introduces in a natural way dark matter components. 
Dark matter in inflationary models which accompanies ordinary matter to govern the evolution of the universe 
can be found in |28[ ]. 

4. We have formulated the conditions for a solution to be that of general relativity, and showed that the Einstein 
case can be achieved only as a limiting case. 

5. Canonical quantization is possible with the aid of Dirac brackets. The resulting Wheeler de-Witt equation 
includes operators which are not free, but are constrained by the second class constraints as operator identities. 

6. The ground is ready for functional integral quantization, the BRST generator is of rank 1, and the determinant 
of second class constraints has been brought to a simple form. 

7. A simple application of geodetic brane gravity to cosmology is possible within the framework of a minisuperspace. 
Classical cosmological models appear in [^5|,^6| . Canonical quantization appears in [^7j , and the complementary 
functional integral quantization in p9| . 

8. Another significant advantage of GBG over GR is the solution to the problem of time. While a homogeneous and 
isotropic metric is characterized by only one dynamical variable (the scale factor of the universe), the embedding 
vector contains two dynamical variables (the scale factor and the bulk time). Thus, taking the embedding vector 
to be the canonical variables, will enhance the theory with one extra variable that may be intepreted as a time 
coordinate. 



SUMMARY 



APPENDIX A: FUNCTIONAL DERIVATIVES 



• Let F[y] be a functional of y(x) such that 8F = J d 3 xf(x)5y(x) then the functional derivative is 



SF 



/(*)• 



5y(x) 



The chain rule holds for functional derivatives — — = — — — 

Sy(x) oG oy(x) 

• The delta distribution is a scalar density of weight 1 such that for a 3-scalar f(x) 





(Al) 



The covariant derivative of the delta function S 3 ^(x — z) is defined for a 3- vector g l (x) 




(A2) 



• The delta function is symmetric with its two arguments 
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8(x — z) — S(z — x) . 

The first covariant derivative of the delta function is antisymmetric with its arguments 

5\i(x - z) = V x ,8(x - z) = -V z i8{z -x) = -8\i(z - a;) . 

While the second covariant derivative is again symmetric. 
• The basic functional derivatives are 



5y A {x) _ , 



Sy B (z) 
8y B {z) 

5y%{x) 

Sy B (z) 



6b5(x - z) 



= (si - yfah ab y B \b) Sftix -z)- y B \ v yt a h ak 5 lk (x - z) 



For a general expression &(x, y, yu,y\ij) the functional derivative is 
Sy^z) = dy^ {x)5{x - Z)+ OyJ {x)S ^ X - Z) 



+ p-W ( § a - yfb hab yA\ a ) s {ij (x -z)- yA\ijy%h hk s [k (x - z) 

V\ij 



Another nontrivial example is the 3-dimensional Christofel symbols T kl = h %d y A yA,ki 



snM 

5y A (z) 



h l3 yA\ki(x)S\j{x -z) + K' :, y A \j{x)8\ k i{x - z) 



• The Poisson brackets are defined in the usual way 

SF 8G 



{F,G} = J d 3 



5F 8G 



8y A (x) 8P A {x) SP a (x) 8y A {x) 



APPENDIX B: POISSON BRACKETS OF CONSTRAINTS 



We will start with the constraints (£> 



o = 



8nG 
2Vh 



8^G 



y\ k -P~o 

Px~0, 



95 = 



8ttG 
2y/h 



8^G 



(A + K {3) ) + P0(tf - Xiy-QP 



P9(* - XI)- 2 QP 



The PB of these constraints are listed below 



{Mx), Mz)} = [Q l (x) + Q l {z)] 8 {i (x - z) « . 
{<f>o(x), cpi(z)} = (/> (z)S\i(x - z) - <j> 5 X,i{z)S(x - z)«0. 
{(f>o(x), (j)i(z)} = (j> 5 (z)8(x -?)«0. 

{Mx),Mz)} = [B\x) + B l {z)] 8 {i (x -z) + M(z)S{x - z) 
{4>k(x),(j>i(z)} = (f>i(x)S\ k (x - z) + 4>k{z)8\i{x - z) s» . 
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{(/>k{x), 




{4>k(x), 




{4>iix), 




{4>i{x), 




{<p5 (x), 





ox' 



dX 



5(x — z) . 



),&(*)} = [F l (x) + F\z)] 6\i(x-z) . 
Where the shorthanded expressions are 



dX 



M 



^£ [pe(* - A/)- 3 ep] 

Vh 



Vh 



Q l = h l3 <f> j +2 
B i 



{Kh ij -K ij ) {j - ^-hvfa 
Vh 

d(j>5 _ 
dX ' 



Vh 

Vh 



'5-0 

dX 



8irG 



d(j) 5 
dX 

X—-K+ (R t3 2K il K l i ) — (Khv - K») 



dX 

+(Kh lJ -K%j 
8irG 



dX 



3'dx 

-2^[P(*-A7)- 1 ] K (* 

\i vh 

d 



_ P(* - A/)- 1 ^ - A/)- 1 ^,,—^ - ir->) 



ia 2 



3 <9A 2 

90s 
dX 

8wG 



(Kh ij - K ij 
dX 



>\j 

-1 



dX 



(Kh ij ~ K ij ) 



d 



-2^P(* - XI)- 2 [(* - A/) _1 P] ! . - ^) 



APPENDIX C: MATTER HAMILTONIANS 

Consider here a few simple matter Lagrangians and Hamiltonians, 
• For a cosmological constant , 

^matter — V 9^-^ ! 

T aP = -2Ag a P . 
The corresponding energy /momentum projections are 



The Hamiltonian is simply 

Ti-matter = —^matter = NVh2K = NVhT 7 

• For a scalar field $(x), 
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The momentum II conjugate to $ is given by 

and the corresponding energy/momentum projections are 



Tun = \ + + V , 



The matter Hamiltonian is 



For a vector field -A^x), 



^matter — — 77— \/~ Ijg^ g 1 "* F^Fxa , 
167T 



The momentum II M conjugate to is given by 

n° = o, 
n* = 



kj 



and the corresponding energy /momentum projections are 



2tt 1 
T nn = — ^'Unj + —Wh kl F lk F 3l , 

T n i — -^h kl HkFu ■ 

The Hamiltonian is 

2tt 1 

H = NVhi—WTLilLj + —Wh kl F ik F ol ) + N l W F, - A n\ = NVhT nn + N l VhT m - A U\ . 

h 167T 

In this case the Hamiltonian picks up another Lagrange multiplier A , and an additional constraint 



(C4a) 
(C4b) 

(C5) 

(C6a) 
(C6b) 

(C7) 

(C8a) 
(C8b) 

(C9a) 
(C9b) 

(ClOa) 
(ClOb) 

(Cll) 
(C12) 



APPENDIX D: THE CENTER OF MASS AND RELATIVE COORDINATES 

We will try to make a canonical transformation to the new system. We will use a global pair Y A (t) 7 Pa(*) to 
describe the total momentum and its conjugate coordinate. And, as relative coordinates we will use the directional 
derivatives zf (x) = y A (x) of the field y A (x) . (This is the analog to a discrete system, where the relative coordinates 
are differences between the coordinates of the various particles involved). 
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The variation of the Action with respect to y A (x) is going to be very similar to the variation with respect to hij , and 
therefore will resemble Einstein's equations. The new set of canonical 'coordinates + fields' Y A , P A , z A (x), ir A (x), 
must obey the canonical PB 



{Y^,P S } =8i, 
{Y A ,^ B (x)} = 0, 
{z A (x),P B } = , 
{z A (x),n B {x)} =6£5i6{x-x). 

We will write the transformation from the old set of fields to the new set as 



Y A (t) = J d 3 xf(x)y A (t,x) 



d 3 xP A (t,x) 



zt(t,x)=y A i (t,x) 



d 3 xP A (t,x)J i (x,x) 



While the inverse transformation is 



y A (t,x) = Y A (t)+ J d 3 xz A (t,x)J l (x,x) 
P A (t,x)=P A (t)f(x)-n A At,x) 



(Dla) 
(Dlb) 
(Die) 

(Did) 

(D2a) 

(D2b) 
(D2c) 
(D2d) 

(D3a) 
(D3b) 



The functions f(x), J l (x,x) are distributions over the V3 manifold, the y do not depend on the canonical fields, and 
in particular are independent of the 3-metric. The solution to Eq.(Dl) put some restrictions on f(x),J l (x,x), and 
they must fulfill the following relations 



d 3 xf{x) = 1 



d 3 x f(x)J l {x,x) = 



dJ l (x, x) 

dx~i 
dJ l (x, x) 

dx % 



= SjS(x - x) , 
= f{x) -S(x-x) 



(D4a) 
(D4b) 
(D4c) 
(D4d) 



We assume one can find such distributions and we move on to the dynamics. We will start with the Hilbert action 
(Pj) and do the variation with respect to the new variables 



SS 



16ttG 

-2 
16ttG 
2 

16ttG 
2 

16 



d A x^(G^ - SnGT^)8g^ 
d 4 x y/=g(G^ - &nGT^)y A>li 5y A 

2 — J d A x^—g(G^ -^GT^)y A ^8z A {x) , 



SY A + / d 3 x8zf(x)J i (x,x) 



The variation with respect to Y A will lead to the conservation of the total momentum 

-2 



16ttG 



d 3 x [V^giG" - SnGT^yA^] — ji A = const. 



(D5) 



(D6) 



The variation with respect to z A {x) will lead to an equation similar to Einstein's equations, but, the right hand side 
does not vanish 
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gy A , a [G m ~ 8irGT at ](x) = / d 3 xJ l (x,x) [V=g(G a0 ~ 87rGT a0 )y A , a (x)] 



Multiply Eq.(|D7|) by ■j=°g' w y* to get 



-9 



G* ~ ^GT^(x) = D<*(x) = —=g^y A v {x) \ d 3 x,P(x,x) [V^(G a0 - ^GT a0 )y A>a {x)] Q 



(D7) 



(D8) 



An Einstein physicist will interpret Eq.( |D8| ) as if there is some additional matter in the universe, and may call it dark 
matter. 

It is easy to reveal Eq.(|^) if one takes the derivative of Eq.(D7) with respect to x % and use (D4d). 



APPENDIX E: DETERMINANT OF SECOND CLASS CONSTRAINTS PB 



We would like to calculate the determinant of C mn (x, z) (}49|). First we will find the eigenvalues of C. Take v(x) to 
be a two components scalar function 



v(x) = 



The eigenvalue equation of C is 



d 3 z C(x, z)v(z) — av(x) 



(El) 



(E2) 



Inserting Eq.(|49|) into Eq.(E2) one can see that the components of v(x) are proportional, and must obey a differential 
equation 



1505, 

g = nTJ ' 



Multiplying ( |E3t| ) by / one gets 



(E3a) 
(E3b) 

(E4) 



Eigenvalues of a differential operator are determined by the boundary co ndit ions. Our boundary conditions are 
actually the fact that the 3-manifold has no boundary. Thus integrating Eq.(E4) over V3 gives us 



d 3 ;, 



1 

a \ o\ 



(E5) 



Arranging Eq.(E5) one gets 



Jd 3 x(^) 2 fHx) 
fd 3 xf 2 {x) 



(E6) 



• G mn (x, z) is a PB matrix and therefore antihermitian, this causes the eigenvalues of C to be purely imaginary. 

• One can see that the eigenvalues of C are affected only by the off diagonal terms 4^-, not by F l . 
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The structure of a is very simple. Define the probability density 

™ - j0h) ' (E7) 

one sees that any eigenvalue of C is simply the expectation value of f with respect to some probability distri- 

bution / 

2 _ (d^ 2 



For each / one finds 2 complex conjugate purely imaginary eigenvalues. The determinant of C is therefore the 
multiplication of all eigenvalues 

d- c -n((^)') / m 

The probability density over a compact manifold can be parameterized by the appropriate harmonics, and the product 
is countable. See for example f30|,pl| for the compact S3 harmonics. 
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